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Executive  Summary 


Mide  Technology  Corporation  was  approached  by  the  U.S.  Army  Research  Laboratory  (ARL)  to 
investigate  the  effect  of  using  the  unique  effects  of  piezoelectric  materials  in  order  to  reduce  the 
magnitude  of  the  ballistic  shock  of  ballistic  events.  The  program  was  set  up  to  develop  a  three- 
dimensional  (3-D)  structural  analysis  code  to  simulate  the  effect  of  the  piezo  materials  on  the 
displacement  of  a  Kevlar1  29  laminate  armor  plate.  These  models  were  to  be  calibrated  via 
simple  ballistic  tests.  Once  the  code  was  calibrated,  an  optimization  and  design  effort  would  be 
completed  in  order  to  perform  proof-of-concept  tests. 

Upon  receiving  the  orthotropic  material  properties  of  the  Kevlar  29  composite,  Mide’s  research 
and  development,  which  specializes  in  smart  materials,  decided  to  use  a  MATLAB2-based  finite 
element  approach  that  accurately  models  the  3-D  anisotropic  behavior  of  the  armor.  The  model 
also  included  the  ability  to  model  piezoelectric  material  and  its  interaction  with  the  Kevlar 
material. 

Some  delays  in  the  process  forced  the  testing  to  be  postponed  until  very  late  in  the  program. 

This  reduced  the  amount  of  time  available  for  the  calibration  and  optimization  of  the  model.  As 
an  alternative,  Mide  agreed  to  integrate  some  of  its  off-the-shelf  packaged  piezoelectric  actuators 
with  Kevlar  composite  targets  in  order  to  test  the  effects  of  the  ballistic  events  on  the  materials. 
This  step  was  scheduled  for  the  Phase  I  Option,  and  so  certain  tasks  from  Phase  I  moved  into  the 
option  and  vice  versa. 

Further  difficulties  in  obtaining  reliable  data  from  the  ballistic  tests  hindered  calibration  of  the 
model.  At  this  stage  it  was  decided  to  simulate  the  tests  that  were  conducted  in  ANSYS3,  a 
commercial  finite  element  analysis  package.  Correlation  between  the  model  and  high  speed  video 
data,  as  well  as  data  from  the  piezoelectric  tests,  was  obtained  as  far  as  possible.  The  model  was 
then  used  to  investigate  the  effect  of  creating  a  destructive  surface  wave  on  reducing  the  maximum 
displacement  of  the  impact  point.  Early,  un-calibrated  and  un-optimized  results  indicate  that  a 
27%  reduction  in  the  displacement  of  the  impact  point  is  possible  when  a  high  voltage  pulse  is 
delivered  to  a  certain  arrangement  of  piezoelectric  materials. 

These  results  are  encouraging  but  should  be  validated  by  testing.  Calibration  of  the  model  is 
scheduled  early  in  the  option.  The  model  is  ready  to  be  optimized  as  soon  as  the  calibration 
parameters  are  updated.  This  optimization  will  assist  in  a  rapid  design  cycle  to  ensure  that 
testing  of  an  integrated  system  will  still  be  performed  during  the  option. 


Kevlar  is  a  registered  trademark  of  E.  I.  DuPont  de  Nemours  &  Co.,  Inc. 
2MATLAB  is  a  registered  trademark  of  The  MathWorks. 

3ANSYS,  which  is  not  an  acronym,  is  a  registered  trademark  of  ANSYS,  Inc. 
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The  following  is  a  summary  of  the  achievements  under  each  of  the  tasks  as  outlined  in  the  work 
breakdown  structure  of  the  Phase  I  proposal. 

Task  1:  Develop  and  Understand  Requirements 

Mide  had  a  productive  “kick-off’  meeting  with  the  technical  monitor  at  Aberdeen  Proving 
Ground  (APG).  The  meeting  restated  the  objectives  of  the  program,  and  the  schedule  and  plan 
for  the  coming  months  were  agreed  upon.  The  technical  monitor  also  provided  Mide  with 
literature  and  defined  test  parameters.  Mide  further  received  detailed  specifications  for  the 
Kevlar  plate  to  be  used  during  the  program. 

Task  2:  Parametric  Modeling 

A  3-D  finite  element  model  was  constructed  to  study  the  anisotropic  behavior  of  the  armor 
material.  This  code  was  verified  with  previously  published  results  and  was  used  to  understand 
the  behavior  of  the  armor  system. 

Task  3:  Absorber  Optimization 

This  task  involved  the  integration  of  piezoelectric  materials  into  the  model.  The  voltage  output 
and  visual  displacement  correlations  from  the  integrated  model  were  validated  with  the  limited 
experimental  data  available.  Actuation  control  techniques,  piezoelectric  material  size  and 
placement,  piezoelectric  shunting  techniques  and  other  shock  mitigation  schemes  need  to  be 
explored  in  a  post-Phase  I  option. 

Task  4:  Preliminary  Experiments 

Ballistic  experiments  using  polypropylene  bullets  shot  at  Kevlar  targets  at  1000  m/s  were 
performed  at  APG.  These  tests  include  shots  at  targets  with  piezoelectric  materials  integrated 
into  them.  Some  useful  data  were  obtained  during  the  tests,  and  more  tests  are  scheduled  in  the 
period  between  the  end  of  Phase  I  and  the  start  of  the  option. 

Task  5:  Shock  Absorber  Design 

A  piezoelectric  material  configuration,  size  and  actuation  scheme  was  chosen  to  investigate  the 
effect  of  producing  a  destructive  surface  wave  in  the  Kevlar  to  reduce  the  displacement  of  the 
impact  point.  The  result  of  this  modeling  effort  was  a  27%  reduction  in  the  translation  of  the 
impact  point.  This  result  is  encouraging,  but  it  must  be  noted  that  the  model  is  not  calibrated  and 
that  the  result  can  only  be  verified  with  experiments  in  a  post-Phase  I  option. 
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1.  Description  of  Work  Completed  During  Phase  I 


This  section  reports  about  the  progress  that  was  made  under  each  of  the  proposed  tasks.  The 
original  tasks  are  listed  in  the  beginning  of  each  chapter  in  italics. 

1.1  Taskl:  Develop  and  Understand  Requirements 

In  this  task,  Mide  worked  with  the  U.S.  Army  Research  Laboratory  (ARL)  to  define  a  set  of  shock 
absorber  requirements.  ARL ’s  experience  and  understanding  of  the  dynamics  of  a  ballistic 
impact  were  used  to  develop  these  requirements  and  to  transfer  knowledge  to  Mide  personnel. 

1.1.1  Kick-off  Meeting 

Mide  had  a  productive  meeting  with  the  technical  point  of  contact  (TPOC)  at  Aberdeen  Proving 
Ground  (APG),  Maryland.  The  meeting  focused  on  the  scheduling  of  tasks  and  the  definition  of 
the  problem. 

1 . 1 . 1 . 1  Problem  Definition 

Current  bulletproof  vest  systems  include  a  ceramic  plate  with  a  woven  Kevlar  29  backing.  The 
ceramic  plate  breaks  up  the  bullet  and  reduces  it  to  small  pieces.  The  Kevlar  backing  then  catches 
these  pieces  of  ceramic  and  bullet,  thus  preventing  them  from  penetrating  the  skin  of  the  Soldier. 

Mide  gained  a  better  understanding  of  the  mechanics  of  the  ballistic  protection  system.  A  main 
concern  for  personnel  protection  is  the  amount  of  displacement  at  the  back  of  the  armor  system 
after  a  ballistic  event.  This  displacement  can  cause  damage  to  organs  if  not  contained,  often 
resulting  in  incapacitation  the  Soldier. 

This  displacement  is  a  direct  result  from  arresting  the  mass  of  target  and  penetrator  that  travels 
through  the  armor  system  after  a  ballistic  event.  The  main  focus  of  this  project  was  to  reduce  the 
amount  of  displacement  at  the  back  of  the  armor  protection  system.  Mide  realized  that  the 
ballistic  event  has  two  components  that  need  to  be  addressed. 

First,  the  energy  of  the  shock  wave  needs  to  be  dissipated.  Piezoceramics  have  been  successfully 
applied  to  a  number  of  high  frequency-damping  problems.  Part  of  the  program  focused  on  how 
to  effectively  employ  these  ceramics  so  that  maximum  energy  dissipation  and  damping  can 
occur. 

The  second  component  is  the  momentum  translation  from  the  high  speed,  lightweight  bullet  to 
the  slowly  moving,  higher  mass  armor  vest.  Ideally,  the  armor  system  translates  all  the  momentum 
into  a  very  slow  moving  mass  that  contacts  a  wide  area  of  the  body,  thereby  reducing  the  pressure 
on  the  body  and  minimizing  the  displacement  of  the  body  tissue. 
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A  successful  armor  protection  system  will  deal  effectively  with  both  these  components  in  order 
to  reduce  the  damage  to  the  body. 

The  TPOC  also  supplied  Mide  with  information  about  previous  ballistic  tests  with  ceramic  and 
composite  targets,  evaluation  of  different  periodic  cellular  materials  and  blunt  trauma  effects 
after  low-velocity  ballistic  impacts. 

1.1.2  Test  Setup  for  Preliminary  Tests 

The  approach  to  Task  1  was  as  follows: 

1.  Non-perforating  small  arms  ballistic  experiments  that  focused  on  collecting  shock  wave 
and  displacement  data  off  the  back  face  of  the  target  were  conducted.  These  preliminary 
tests  confirm  the  non-perforating  test  velocity  for  the  test  projectile  into  the  target.  The 
data  that  were  collected  are  as  follow: 

•  Back  face  displacement  versus  time  for  0-  to  150-kHz  range. 

•  Shock  wave  speed  and  possibly  shock  wave  pressure. 

•  Maximum  back  face  displacement. 

2.  Another  experimental  setup  used  was  the  composite  panel  attached  to  a  steel  plate  and 
measured  the  force  transmitted  during  the  projectile  impact  with  a  load  cell.  Here  would  be 
an  example  of  a  worst  case  scenario  of  the  force  transmitted  to  the  panel  because  of  the 
coupling  of  the  steel  plate. 

3.  The  data  were  obtained  with  the  following  transducers: 

•  A  crush  gauge  (or  piezo  pin  transducers)  with  signal  conditioning  was  initially  used  to  get 
an  initial  measurement  of  the  maximum  back  face  displacement  off  the  back  of  the  target. 

•  An  eddy  current  gauge  was  set  up  at  the  back  of  the  panel  to  give  a  displacement  versus 
time  reading  and  works  for  the  0-  to  150-kHz  frequency  range.  Velocity  was  obtained 
through  post-processing  differentiation. 

•  A  polyvinylidene  fluoride  transducer  was  used  to  measure  the  shock  wave  speed  from  the 
front  face  to  the  back  face  of  the  target  and  the  shock  wave  pressure. 

4.  The  data  acquisition  system  was  four,  four-channel  digital  oscilloscopes. 

1.1.3  Preliminary  Test  Definitions 

The  properties  of  the  materials  that  were  used  in  the  preliminary  tests  and  some  of  the  test 
parameters  are  listed  in  table  1 . 
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Table  1.  Material  properties  of  the  armor  system. 


Parameters 

Target  -  Kevlar 

Units 

Value 

Volume 

cm3 

331.84 

Size  (L  ,  W,  t) 

cm 

mm 

L=W=22.86 

t=6.35 

Material  Density 

g/cm3 

1.44 

Poisson’s  Ratio  (in-plane) 

-- 

0.12 

Poisson’s  Ratio  (out-of-plane) 

-- 

0.4 

Material  Stiffness  -  x  and  y  directions 

GPa 

24.1 

Material  Stiffness  -  z  direction 

GPa 

10.4 

Fiber  Shear  Strength 

GPa 

0.24 

Matrix  (Delamination  Shear)  Strength 

MPa 

48 

Material  Compression  Strength 

GPa 

0.35 

Thru  Thickness  Crush  Strength 

MPa 

69 

Penetrator  -  Polypropylene  Plastic 

Right  cylinder  projectile 

Units 

Value 

Diameter 

mm 

12.57 

Length 

g 

14.78 

Weight 

g 

1.56 

Density 

g/cm3 

0.91 

Velocity  (at  impact) 

m/s 

300-up 

Tensile  strength  at  yield 

MPa 

35.85 

Tensile  elongation  at  break 

% 

600 

Tensile  modulus 

GPa 

1.65 

Flexural  modulus 

GPa 

1.59 

Flexural  strength 

MPa 

48.26 

Notched  izod  impact  strength 

J/m 

60 

1.2  Task  2:  Parametric  Modeling 

In  this  task,  Mide  worked  with  ARL  to  develop  a  parametric  model  with  which  the  performance 
of  a  shock  absorber  design  was  evaluated  and  validated. 

Upon  receiving  the  orthotropic  material  properties  of  the  Kevlar,  Mide  decided  to  use  a  MATLAB4- 
based  finite  element  approach  that  accurately  models  the  3-D  anisotropic  behavior  of  the  armor. 

The  code  is  attached  as  appendix  A.  This  code  was  verified  as  far  as  possible  and  it  was  used  to 
understand  the  dynamics  of  the  impact  and  the  resultant  stresses  that  a  piezo  mitigation  system 
must  address.  Figure  1  shows  the  model  and  figures  2  through  4  are  the  first  three  flexible  modes 
predicted  by  the  model.  Results  of  an  impact  simulation  are  in  figures  5  through  7  and  time 
responses  in  figures  8  and  9. 


4MATLAB  is  a  registered  trademark  of  The  MathWorks. 
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Mode  Number:  7,  Frequency:  306.136  Hz 


Figure  2.  First  non-rigid  resonant  mode  shape. 
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Mode  Number:  9,  Frequency:  543.986  Hz 


Y-Direction  (m) 


X-Direction  (m) 


Figure  4.  Third  non-rigid  resonant  mode  shape. 
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Figure  5.  Von  Mises  stress  at  0.8  |lis  after  the  1.56-gram  projectile  impacted  the  center  of  the 

plate  at  450  m/s.  (Note  that  the  top  surface  is  not  deformed  yet  and  that  the  damage  is 
limited  to  the  impact  point.) 


100 

— i — 

i 


150  200  250  300 


Figure  6.  Von  Mises  stress  at  4  ps  after  the  1.56-gram  projectile  impacted  the  center  of  the  plate 
at  450  m/s.  (The  stress  on  the  top  surface  has  increased  but  little  or  no  deformation  is 
visible.  The  maximum  top  surface  Von  Mises  Stress  is  around  155  MPa.) 
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VM  Stress  in  MPa  -  Step:  200/4000,  Time=  8.0  micro-sec.,  Scaled  by=1 .0 


50  100  150  200  250  300  350 


Figure  7.  Von  Mises  stress  at  8  |Lis  after  the  1.56-gram  projectile  impacted  the  center  of  the  plate 

at  450  m/s.  (Clear  deformation  is  visible  at  this  point  with  the  top  surface  near  350  MPa.) 
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Figure  8.  Displacement  and  velocity  (deformation)  of  impact  point. 
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Figure  9.  Displacement  and  velocity  (deformation)  of  top  surface  directly  above 
the  impact  point. 


1.3  Task  3:  Absorber  Optimization 

The  models  of  Task  2  were  used  to  establish  a  design,  to  optimize  the  design,  and  to  estimate  the 
performance  of  the  design. 

1.3.1  Finite  Element  Model  Formulation 

Mide  uses  the  National  Center  for  Supercomputing  Applications  structural  analysis  system 
Nastran5,  which  is  a  finite  element  analysis  program.  Nastran  does  not  have  a  piezoelectric 
element  in  its  library.  Initially,  the  cost  of  ANSYS  which  does  have  a  piezoelectric  element  was 
considered  prohibitive,  so  Mide  decided  to  evaluate  the  performance  of  different  geometries  of 
active  (piezo)  and  passive  materials,  with  its  MATLAB-based  finite  element  code.  It  is  easier  to 
generate  many  models  in  MATLAB  compared  to  the  alternative  of  using  Nastran  to  generate  the 
finite  element  models  and  then  importing  the  finite  element  models  into  MATLAB  where  the 
active  elements  can  be  modeled. 

The  Mide  MATLAB  finite  element  code  is  based  on  the  hybrid  stress  finite  element  that  has 
excellent  accuracy  even  if  the  elements  are  distorted.  This  sub-section  provides  a  quick  review 
of  the  formulation  of  hybrid  stress  finite  elements. 


5Nastran  is  a  trademark  of  the  National  Center  for  Supercomputing  Applications. 
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In  the  hybrid  finite  element  model,  the  strain-stress  relationship  is 

£  -  S<J 

The  stress  is  expressed  in  terms  of  the  undetermined  stress  parameters  p : 

a  =  pp 


From  the  strain-displacement  ( q )  relations, 

s  -  Bq 

The  Hellinger-Reissner  functional  ( ff ) 

^-R=\\rB{(7‘i)+(T‘Js‘j]dv~\uiTids 

v  st 


With  nR=~pTHp  +  pTGq-(Bq)TQ 

where  Q  are  the  external  forces  and  where 


The  operations  aV\R  /80r=  0  yield 

Eliminating  p  and  substituting 

we  obtain 

where  the  stiffness  matrix  is  given  by 


H  =  J  PT  SPdV 

V 

G  =  J  BTPdV 

V 


P  =  H'Ga 


a  =  B  'q 


kq  =  Q 


where 


h  =  r  pt  s  p  i  w 

11 18x18  18x6^6x6^  6x18  J  nl 


^  —  ^24x18^18x18^18x24  : 


Intg 


B  P 

**  24x6^6x18 


-i  PT  R 
VpT  o  p  1  _J_8x6_j 

|_"M8x6^ox6  cixl8  J  j 


J 


6x24  W 

nintg 


GnxU=\BTPdV  = 


P  B 

1  18x6^6; 

kl 


-W, 


Intg 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


(11) 


(12) 


1.3.2  Anisotropic  Materials 

The  first  modification  of  the  existing  MATLAB  code  was  to  include  in  the  model  the  ability  to 
use  anisotropic  materials.  Given  the  Jacobian  of  an  element, 
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The  orientation  of  element  in  the  global  coordinate  system,  given  its  nodal  connectivity  numbering, 
can  be  determined  as  the  normalized  inverse  of  the  Jacobian: 


Given  the  definition  of  the  compliance  matrix  S , 
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The  compliance  matrix  S  in  the  element  coordinate  system  is 
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1.3.3  Modeling  of  Piezoelectric  Materials 


The  next  step  in  expanding  the  capabilities  of  the  MATLAB  code  was  to  determine  the  voltage 
generated  by  the  piezo  elements  when  deformed  and  the  deformation  they  induce  when  a  voltage 
is  applied.  A  piezoelectric  material  is  transversely  isotropic.  In  the  Institute  of  Electrical  and 
Electronics  Engineers  Standard  (STD- 176- 1978),  the  relationship  among  dielectric  displacement 
(D),  mechanical  stress  (T),  the  electric  field  (E)  and  mechanical  strain  is  (S) 
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where  s  is  the  permittivity  of  the  material  and  it  is  assumed  to  be  isothermal.  In  engineering 
notation, 
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This  can  be  re-written  in  block  matrix  form  as 
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where  the  coupling  coefficient  is  now  e  and  c  represents  stiffness. 
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The  piezo  can  be  a  sensor  (2  =  0),  which  yields  the  following  equations: 

/t  =  sE  +  ecj  =  0  =>  E  =  -(f)  *  ecr 
s  =  eT E  +  SE<j  =  -eT  (£)  1  ea  +  SEa 

or  it  can  be  driven  by  an  applied  electric  field: 

^  =  ZE applied  +e0’ 


(22) 


(23) 


^ E appi^d  is  a  forcing  term  and  can  be  included  in  the  finite  element  model  as  n  virtual  work  terms. 
This  was  implemented  in  the  MATLAB  code. 


Place  this  in  the  Hellinger-Reissner  functional: 


ff  =  p  +  j  <JveTEApplieddV  +  /3tJ  a  -  (Baf  Q 

1  V 

1 

(24) 

=  -^j3THj3  +  j3TF  +  /3TJa-(Ba)T  QJ3T 

Where  F 

F  E  =  f  PT  eT  dVE 

(25) 

This  yields  a  new  “forcing”  term: 

QI./J'F£  =  [B-,Lp-lL,: 
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The  active  strain  is 


(27) 


2  —  ^ixi^ix\  +  eix(P(>x\ 


when  there  is  an  applied  field  and  without  an  applied  field,  the  generated  field  is 

A  =  £E  +  ecr  =  0=>E  =  -(£ )  '  ecr 
£  =  eTE  +  SEcr  =  -eT  (s)-'  ecr  +  SEcr 


(28) 


Note  that  the  same  coordinate  transformation  needs  to  be  performed  if  the  material  properties  are 
aligned  along  the  element  coordinate  system. 

Figure  10  is  a  simple  example  where  a  voltage  was  placed  over  an  element  in  the  Z-direction  (the 
<7,,  direction). 


Figure  10.  Deformation  of  a  piezo  cube  when  a  voltage  is  applied 
across  the  c/33  (Z-)  direction.  (The  bottom  comers  are 

fixed  and  the  cube  expands  in  the  Z-direction  while 
contracting  in  the  X-  and  Y-directions.  The  constraints 
introduce  stresses  in  the  cube.) 


1.3.4  Linear  Results 

The  finite  elements  were  programmed  in  FORTRAN  (Formula  Translator)  and  linked  with 
MATLAB  executable  files.  The  model  was  used  to  examine  the  performance  of  the  ceramic-glass- 
reinforced  plastic  (GRP)  composite  targets  of  Strapburger,  Senf,  Burkins,  and  Gooch  (1996).  This 
section  reports  the  results  where  the  material  behavior  was  assumed  to  be  linear.  The  next  step  was 
to  perform  a  nonlinear  large  deformation  analysis  where  elements  in  which  the  stresses  exceed 
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failure  criteria  are  eliminated  from  the  analysis.  Figure  1 1  shows  the  linear  first  mode  predicted  by 
the  Mide  code.  Note  that  for  these  results,  clamped  boundary  conditions  were  assumed  in  the 
model.  The  total  target  mass  is  predicted  to  be  1 1.40  kg  or  25.08  lb. 

The  model  was  also  used  to  examine  the  behavior  of  the  target  if  the  boundary  conditions  were 
assumed  to  be  simply  supported. 


Figure  11.  First  mode  shape  of  the  ceramic-GRP  composite  target  predicted  by  the  Mide  code. 
(Boundary  conditions  were  clamped.) 
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Figure  12.  Displacement  and  velocity  (deformation)  of  impact  point.  (Boundary 
conditions  were  clamped.) 
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Figure  13.  Displacement  and  velocity  (deformation)  of  the  top  surface  directly  above  the 
impact  point.  (Boundary  conditions  were  clamped.) 


Figure  14.  Predicted  von  Mises  stresses  at  230  microseconds  (ps)  after  impact.  (Boundary  conditions 
were  clamped.) 
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Mode  Number:  1,  Frequency:  383.317  Hz 
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Figure  15.  First  mode  shape  of  the  ceramic-GRP  composite  target  predicted  by  the  Mide  code.  (Boundary 
conditions  were  simply  supported.) 
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Figure  16.  Displacement  and  velocity  (deformation)  of  impact  point.  (Boundary  conditions 
were  simply  supported.) 
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Figure  17.  Displacement  and  velocity  (deformation)  of  the  top  surface  directly  above  the 
impact  point.  (Boundary  conditions  were  simply  supported.) 


Figure  18.  Predicted  von  Mises  stresses  at  230  micro-seconds  after  impact.  (Boundary  conditions 
were  simply  supported.) 
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As  expected,  the  linear  analysis  was  substantially  different  from  the  results  reported  in  Strapburger 
et  al.  (1996).  However  Mide  code  can  now  predict  the  voltages  generated  by  piezo-materials  and 
is  ready  for  studies  of  how  the  piezo-elements  can  be  used  to  reduce  the  deformation  of  the  back 
plate. 

1.3.5  ANSYS  Finite  Element  Model 

Initially,  Mide  considered  the  cost  of  obtaining  ANSYS  specifically  for  this  program  to  be 
prohibitive.  Another  Mide  program  had  the  funds  available  and  decided  to  purchase  the  FEA 
package.  This  allowed  Mide  to  generate  an  alternate  finite  element  model  for  this  program.  The 
advantages  of  using  ANSYS  instead  of  the  MATLAB  model  are  as  follow: 

•  A  finer  grid  can  be  modeled. 

•  The  actual  models  can  be  created  in  a  solid  modeling  program  and  imported  into  the  FEA 
package. 

•  Optimization  is  supported. 

•  Highly  detailed  piezoelectric  models  are  available. 

•  Memory  management  schemes  in  ANSYS  allow  for  more  detailed  analysis. 

•  Time  history  and  animation  plots  are  easily  obtained. 

It  was  thus  decided  at  a  late  stage  in  the  program  to  change  to  the  ANSYS  model. 

1.3. 5.1  ANSYS  Model  Description 

Figure  19  illustrates  the  meshed  ANSYS  model.  The  model  consists  of  tetrahedral  elements  for 
the  Kevlar  parts  and  solid  bricks  for  the  piezo  parts.  Figure  19  shows  the  model  in  a  deformed 
state.  The  model’s  boundary  conditions  are  placed  at  three  nodes  on  every  comer  on  the  piezo 
side  of  the  Kevlar  plate  in  every  direction.  A  simple  assumption  was  to  model  the  bullet  impact 
as  an  initial  velocity  on  the  nodes  surrounding  the  center  of  the  plate  with  a  0.5-inch  diameter. 

This  velocity  was  set  at  1000  m/s. 

The  piezo  elements  were  constrained  so  that  all  the  elements  on  their  top  and  bottom  (attached  to 
the  Kevlar)  surfaces  were  electrically  connected  to  each  other.  This  had  the  same  effect  as  placing 
an  electrode  over  these  surfaces  and  modeled  the  QuickPack6  and  PowerAct  sensors  and  actuators 
well. 

A  typical  result  of  the  analysis  is  shown  in  figure  20.  Here,  the  elements  are  deformed  and  high¬ 
lighted  according  to  their  displacement  in  the  z  direction  (direction  the  bullet  is  traveling).  The 
shock  waves  can  clearly  be  seen  on  the  surface,  causing  the  displacements. 


6QuickPack  and  PowerAct  are  trademarks  of  Mide  Technology  Corporation. 
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Figure  19.  The  meshed  ANSYS  model. 
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Figure  20.  A  typical  result  from  the  ANSYS  analysis  illustrating  the  nodal  displacement 
in  the  z  direction. 
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1.3. 5.2  Time  Step  Displacement  Comparisons 

The  following  series  of  images  compares  the  displacement  output  of  the  AN SYS  model  to  that  of 
the  high  speed  video  camera  for  test  2649.  The  high  speed  video  camera  has  a  shutter  speed  of 
121  ps/frame.  The  model  was  run  with  a  10- ps  interval.  Every  12th  step  of  the  output  of  the 
model  is  shown  in  table  2.  The  plate  is  c-clamped  to  the  holder  and  struck  on  the  face  contacting 
the  frame.  The  impact  is  below  the  black  tape. 

It  is  interesting  to  note  how  the  model  captures  the  initial  ballistic  event,  causing  the  sharp  bubble 
in  the  Kevlar,  and  the  secondary  wave  that  displaces  the  Kevlar  from  the  frame.  This  comparison 
shows  that  the  model  is  somewhat  representative  of  reality.  Unfortunately,  real  displacement  data 
were  not  available  and  a  true  calibration  of  the  model  was  not  possible. 


Figure  21.  Comparison  of  high  speed  video  and  the  ANSYS  model. 
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Figure  21.  (continued). 
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Figure  21.  (continued). 
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1.3. 5. 3  Voltage  Generation  Comparison 

This  section  presents  a  brief  comparison  between  the  results  obtained  from  the  only  good  sensor 
(PA1)  in  the  tests  described  in  section  1.4.2.  This  sensor  was  the  only  one  that  produced  data 
with  some  value.  A  part  of  the  sensor’s  output  is  shown  in  figure  22,  where  the  figure  zooms  in 
on  the  first  millisecond  of  data  for  comparison  with  the  ANSYS  model.  The  result  of  the  voltage 
output  of  the  piezo  wafer,  placed  at  the  same  spot  as  the  PA1  sensor  is  shown  in  figure  23. 

Although  the  two  plots  differ  in  many  ways,  there  are  some  remarkable  resemblances: 

•  Both  have  a  dominant  frequency  of  ~6.5  kHz. 

•  The  maximum  voltage  for  the  first  peak  is  ~900  V. 

•  A  similar  high  frequency  content  is  superimposed  on  the  dominant  frequency. 

•  There  is  a  general  increase  in  voltage  over  the  first  1  ms  of  the  event. 

Although  these  similarities  are  vague,  it  is  encouraging  to  realize  that  these  similarities  exist 
without  the  model  being  calibrated.  It  may  be  safe  to  assume  that  trends  that  are  identified  by 
the  model  might  be  achievable  in  reality. 


Piezo  Gauge 


Figure  22.  Voltage  output  from  the  PA1  piezo  sensor  during  the  ballistic  test. 
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Figure  23.  Voltage  output  from  the  PA1  piezo  sensor  in  the  ANSYS  simulation. 


1.4  Task  4:  Preliminary  Experiments 

Small  experiments  were  designed  and  conducted  to  verify  critical  elements  of  the  design.  The 
experimental  results  were  used  to  update  and  verify  the  parametric  models. 

Mide  and  the  ARL  decided  to  combine  testing  of  the  baseline  Kevlar  sheets  and  preliminary 
piezoelectric  sensor  devices.  The  piezoelectric  tests  were  scheduled  for  the  option,  but  since  the 
test  facilities  were  available,  it  was  decided  to  evaluate  the  piezoelectric  devices  during  a 
ballistic  event  and  move  the  design  and  optimization  research  to  a  next  phase  effort. 

1.4.1  Kevlar  Plate  Tests 

Six  tests  were  conducted  to  determine  the  displacement  of  the  Kevlar  sheets  in  response  to  the 
impact  of  a  high  velocity  polypropylene  bullet. 

1.4. 1.1  Test  Setup 

The  test  setup  schematic  is  shown  in  figure  24.  The  controller  initiates  the  event  by  sending  an 
electric  pulse  to  the  gun.  This  pulse  activates  the  firing  pin  of  the  gun  and  ignites  the  charge. 

The  charge  propels  the  polypropylene  bullet  through  the  gun’s  barrel  (figure  26)  and  toward  the 
Kevlar  target  (figure  25).  On  the  way  to  the  target,  the  bullet  passes  through  the  trigger.  The 
trigger  sends  a  signal  to  the  controller  that  commands  the  oscilloscope  and  the  laptop  (high-speed 
camera)  to  start  recording.  The  eddy  displacement  gauge  (figure  26)  records  the  displacement  of 
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the  metal  strip,  and  the  laptop  computer  records  the  sequence  of  events  that  follows  the  impact  of 
the  bullet. 


Figure  24.  Kevlar  displacement  test  setup. 


Figure  25.  Test  setup  as  viewed  from  the  gun  barrel. 
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Figure  26.  The  gun,  left  and  side  view  of  the  eddy  gauge,  metal  strip  and  Kevlar 
target,  right. 


1.4. 1.2  Test  Sequence 

Five  tests  were  performed  on  the  bare  Kevlar  sheets.  Unfortunately,  the  metal  strips  used  were 
not  able  to  withstand  the  ballistic  event.  Only  limited  data  from  these  tests  were  available,  but 
none  of  them  were  useful  to  calibrate  the  model.  Additional  tests  were  mostly  conducted  in  an 
attempt  to  improve  the  performance  of  the  metal  strips  so  that  displacement  data  could  be 
recorded. 

Table  2  summarizes  the  tests  that  were  performed  at  APG.  These  tests  can  be  divided  into  the 
displacement  and  piezo  tests.  The  displacement  tests  were  the  originally  scheduled  tests  for  this 
part  of  the  program.  The  results  from  these  tests  were  to  be  used  for  the  calibration  of  the  FEA 
model.  As  stated  before,  the  tests  ended  mostly  in  the  destruction  of  the  metal  pieces  needed  by 
the  measuring  equipment,  and  thus,  no  usable  data  were  obtained  from  them. 

Figures  27  to  32  show  the  failures  of  the  different  tests.  Initially,  the  metal  strips  for  the  eddy 
current  gauge  were  bonded  to  a  protective  layer  that  covered  the  Kevlar  plate.  The  existence  of 
this  layer  was  only  discovered  after  test  2647.  Although  this  layer  was  removed  for  the  subse¬ 
quent  tests,  the  bonding  layer  between  the  metal  and  the  Kevlar  was  not  strong  enough  to  with¬ 
stand  the  ballistic  event.  Ensuing  tests  will  have  the  metal  integrated  with  the  Kevlar  and  this 
should  solve  the  problem.  The  QuickPack  and  PowerAct  tests  produced  some  useful  data  and 
these  tests  are  described  in  some  more  detail  next. 
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Table  2.  Test  sequence  and  failure  mechanisms  for  a  ballistic  test  conducted  at  ARL. 


Test  No. 

Description 

Failure  Mechanism 

2646 

Displacement  test.  Three  2-inch  metal  tape  strips 

Metal  tape  ruptures 

2647 

Displacement  test.  Metal  sheet  bonded  to  structure 

Metal  sheet  ruptures 

2648 

Displacement  test;  0.75-inch-diameter  metal  tape 
placed  under  the  impact  point. 

Metal  disc  intact;  data  not  useful  for  calibration 
of  model. 

2649 

Displacement  test;  0.75-inch-diameter  metal  disc 
bonded  to  structure. 

Disk  blown  off  structure. 

2650 

Displacement  test;  0.75-inch-diameter  thin  metal 
disc  bonded  to  structure 

Thin  disk  blown  off  structure. 

2651 

Multi  layer  piezo  test. 

Piezos  bonded  to  protective  layer.  Piezos 
delaminated  from  structure. 

2652 

QuickPack  and  PowerAct  piezo  wafer  test. 

Sensors  bonded  to  protective  layer.  QuickPacks 
failed  structurally.  One  PowerAct  intact. 

Figure  27.  Setup  before  test  number  2646. 


Figure  28.  Failure  during  test  2646.  (Note 

how  the  tape  tears  in  all  directions, 
indicating  severe  strain-related 
failure  of  the  metal.) 
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Figure  29.  Failure  of  test  2647.  (Again,  the  rupture  of 
the  metal  sheet  indicated  that  the  strain 
limit  of  the  material  was  exceeded 
significantly.) 


Figure  30.  The  impact  point  is  well  above  the  position 
of  the  0.75-inch  metal  disk  in  test  2648. 
(The  impact  does  not  dislodge  the  tape, 
indicating  that  a  smaller  area  might  have  a 
chance  of  surviving  the  high  differential 
strain  environment  of  the  ballistic  test.) 
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Figure  3 1 .  Delaminating  of  the  0.75-inch  metal  disk 
during  test  2649. 


Figure  32.  Failure  of  test  2650  was  attributable  to  the 
dislodging  of  the  metal  disk. 


1.4.2  Integrated  Piezo  Tests 

Although  tests  with  piezoelectric  elements  are  part  of  the  option  for  this  program,  Mide  decided  to 
combine  these  tests  with  the  Kevlar  plate  tests.  The  main  reason  for  this  is  to  fail  early  and  learn 
lessons  in  order  to  succeed  later.  In  order  to  accomplish  this,  two  sets  of  piezoelectric  sensors  and 
actuators  were  bonded  to  the  Kevlar  surface.  As  stated  earlier,  the  surface  was  covered  with  a 
“peel-off’  protective  layer,  the  existence  of  which  was  only  discovered  during  the  tests  at  ARL. 
Thus,  the  piezoelectric  materials  were  actually  bonded  to  the  protective  peel-off  layer. 
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Two  Kevlar  plates  were  instrumented  with  piezos  for  tests  265 1  and  2652.  The  first  set  of  piezo¬ 
electric  elements  consisted  of  four  layers  of  piezoelectric  arrays.  Each  layer  consisted  of  six  rows 
of  piezos  that  were  made  from  three  separate  1.8-inch  by  0.83-inch  by  0.006-inch  piezoelectric 
wafers,  coupled  in  parallel  to  each  other.  The  second  set  of  piezos  consisted  of  two  sets  of  Mide’s 
PA  ION  PowerAct  and  QP10N  QuickPack  actuators.  The  description  and  summary  of  the  tests  are 
given  in  table  3. 


Table  3.  Description  of  the  integrated  piezo  tests. 


Test  No. 

Description 

Failure  Mechanism 

2651 

Multi  layer  piezo  test. 

Piezos  bonded  to  protective  layer.  Piezos  delaminated  from  structure. 

2652 

QuickPack  and  PowerAct 
piezo  wafer  test. 

Sensors  bonded  to  protective  layer.  QuickPacks  failed  structurally. 

One  PowerAct  intact. 

Because  of  the  presence  of  the  protective  layer,  the  initial  shock  was  enough  to  remove  the  piezos 
from  the  Kevlar,  and  no  useful  data  were  obtained.  It  is  not  deemed  necessary  to  describe  the  test 
and  the  configurations  or  the  results  at  any  length. 


Figure  33.  Four  layers  of  piezoelectric  elements  with  various  electric  configurations  were  bonded  to  the  target 
as  shown  on  the  left.  (Unfortunately,  because  of  the  presence  of  a  protective  layer,  the  piezos  were 
not  bonded  to  the  structure  and  the  failure  [right]  resulted  during  test  265 1 .  No  useful  data  were 
obtained.) 

The  second  piezo  tests  had  four  encapsulated  piezoelectric  wafers  bonded  to  the  Kevlar  plate. 
The  test  setup  is  shown  in  figure  34.  Two  QuickPack  QPlONs  and  PowerAct  PAlONs  were 
placed  opposite  each  other  as  shown.  The  details  of  the  sensors  are  listed  in  table  4.  A  short 
description  of  the  two  types  of  actuators  is  available  at  http://www.mide.com. 
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Figure  34.  Pre-  and  post-test  2652  pictures  of  the  QuickPack  and  PowerAct  sensors. 


Table  4.  Details  of  the  four  encapsulated  piezoelectric  wafers  bonded  to  the  Kevlar  plate. 


Sensor  No. 

Description 

Wafer  Size  (in.) 

Serial  No. 

Piezo  Material 

Capacitance  (nF) 

PA1 

Conformable  wafer 

1.83  x. 83  x0.01 

PA11201 

PZT5A 

41.7 

PA2 

Conformable  wafer 

1.83  x. 83  x0.01 

PA11010 

PZT5A 

45 

QP1 

Rigid  wafer 

1.83  x. 83  x0.01 

QP11205 

PZT5A 

54.4 

QP2 

Rigid  wafer 

1.83  x. 83  x0.01 

QP11214 

PZT5A 

52.6 

Although  the  wafers  were  bonded  to  the  protective  layer,  some  useful  data  were  obtained  from  the 
sensors.  These  results  are  shown  in  figures  35  through  39  and  were  used  in  the  modeling  section 
to  compare  with  the  predicted  output  from  the  sensors.  The  picture  on  the  left  of  figure  34  shows 
the  effect  that  the  ballistic  event  had  on  the  sensors.  PA1  is  the  only  sensor  that  remained  intact. 
All  the  other  sensors  were  destroyed  by  the  event.  It  was  expected  that  the  high  strains  would 
crack  the  rigid  QuickPacks.  The  failure  of  the  one  PowerAct  can  only  be  explained  by  the  fact 
that  it  was  bonded  to  the  protective  layer,  which  delaminated  during  the  event,  causing  the 
PowerAct  to  be  accelerated  away  from  the  surface  and  to  experience  abnormal  bending. 

Figure  35  shows  the  results  from  the  eddy  current  gauge.  The  0.75-inch  metal  strip  was  removed 
form  the  Kevlar  plate  by  the  ballistic  event.  This  caused  a  spike  in  the  gauge  data  and  indicated 
the  time  of  impact.  No  further  useful  data  can  be  derived  from  the  signal. 
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Figure  35.  Results  from  the  eddy  current  gauge.  (This  result  is  indica¬ 
tive  of  the  failure  of  the  metal  strip  to  remain  on  the  Kevlar 
composite  plate.  However,  it  indicated  the  initial  event  at 
0.1  sec  after  trigger,  which  is  useful  for  analysis  of  the 
piezoelectric  signals.) 

Figure  36  shows  the  signal  obtained  from  PA1.  This  is  the  only  useful  piece  of  data  available 
from  the  entire  test  series.  This  trace  was  further  analyzed  in  the  modeling  section  where  it  was 
compared  with  the  ANSYS  model.  Here,  we  can  note  a  few  features  of  the  signal.  The  first 
resonant  mode  is  at  140  Hz.  The  initial  impact  is  sensed  but  does  not  drive  the  signal  to  the 
limit.  The  shock  arriving  around  0.01  second  after  impact  drives  the  signal  out  of  range, 
indicating  that  this  is  the  main  event. 


Figure  36.  Output  from  PA1.  (The  continuous  trace  after  the  event 
indicates  that  these  are  useful  data.  The  spike  at  0. 1 1  sec 
exceeded  the  limits  of  the  oscilloscope.) 
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Figure  37  shows  the  result  form  QP1.  The  impact  cracked  the  piezo  into  a  number  of  smaller 
pieces.  The  data  show  this  as  the  different,  disconnected  pieces  generate  their  individual  signals 
that  are  combined  in  the  trace.  This  combined  signal  appears  as  a  random  trace,  as  can  be  seen 
in  figure  37. 


Piezo  Gauge 

300  | - 1 - 1 - 


Time  (sec.) 

Figure  37.  Output  from  QP1.  (It  is  clear  that  the  sensor  produced  random 
signal  because  the  wafer  was  cracked.) 

Figures  38  and  39  show  what  happens  when  the  wafer  becomes  dislodged  from  the  Kevlar.  The 
ballistic  wave  pushes  the  piezo,  having  momentum,  back.  This  causes  the  element  to  bend  and 
completely  break  off,  as  can  be  seen  in  the  left  picture  of  figure  34.  The  output  voltage  traces 
from  the  elements  indicate  an  initial  charge  that  then  bleeds  off  after  the  piezo  is  destroyed. 


Piezo  Gauge 


Time  (sec.) 


Figure  38.  Output  from  PA2.  (The  sensor  breaks  immediately  after 
the  initiating  event,  bleeding  its  charge  to  the  oscilloscope.) 
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Piezo  Gauge 


Time  (sec.) 


Figure  39.  Output  from  QP2.  (This  sensor  also  has  a  catastrophic  failure, 
thus  the  bleed-off  of  the  initial  signal.) 


1.5  Task  5:  Shock  Absorber  Design 

With  the  parametric  model  developed  in  Task  2  and  updated  in  Task  4,  design  a  detailed 
sandwiched  piezoelectric  shock  absorption  system  that  can  be  fabricated  and  tested  in  the 
option. 

The  ANSYS  model  was  used  to  determine  the  effect  of  sending  a  high  voltage  pulse  to  piezo 
wafers  that  were  placed  around  the  impact  point  in  order  to  create  a  disturbance  wave.  It  was 
hypothesized  that  this  disturbance  wave  might  reduce  the  displacement  of  the  cone  at  the  impact 
point  and  thus  diminish  injuries  and  damage  during  a  ballistic  event. 

The  model  was  used  to  simulate  the  effect  that  0.5 -inch-thick  piezo  wafers  might  have  on  the 
displacement  of  the  Kevlar  composite  plate.  It  is  common  practice  to  limit  the  electric  field  in  a 
piezo  to  20  V/mil.  This  was  the  maximum  field  that  was  used  in  the  model.  Again,  the  initial 
conditions  of  the  model  were  to  constrain  the  comers  of  the  Kevlar  plate  to  zero  movement  and 
to  add  an  initial  velocity  of  1000  m/s  to  the  node  in  a  0.5-inch  diameter  around  the  center  of  the 
plate.  Two  simulations  were  ran,  one  with  the  voltage  pulse  and  one  without  for  comparison. 

Figure  40  shows  the  displacement  of  the  center  node  only  in  response  to  the  ballistic  event.  The 
node  is  displaced  a  maximum  of  230  micro-strains  0.38  ms  after  bullet  impact.  Figure  41  shows 
the  displacement  of  the  same  node,  but  this  time  with  a  destructive  wave  generated  by  the  piezos. 
The  maximum  displacement  is  180  micro-strains  0.39  ms  after  the  ballistic  event.  This  is  a  27% 
reduction  in  displacement.  It  is  further  interesting  to  note  the  high  frequency  element  present  in 
the  piezo-actuated  simulation.  This  high  frequency  element  can  be  responsible  for  the  reduction 
in  the  center  node  amplitude. 
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Figure  40.  Displacement  of  the  center  node  in  response  to  the  ballistic  event. 


Figure  41 .  Displacement  of  the  center  node  in  response  to  the  ballistic  event  when  a 
structural  interference  wave  is  generated  by  the  piezoelectric  elements  in 
response  to  a  high  voltage  pulse. 
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These  are  encouraging  results,  but  caution  must  be  exercised.  The  actual  meaning  of  the  displace¬ 
ment  units  can  only  be  determined  when  the  model  is  calibrated.  This  calibration,  and  only  this 
calibration,  will  quantify  these  results. 

Once  the  model  is  calibrated,  the  AN SYS  model  can  be  used  to  find  the  best  way  to  reduce  the 
displacement  of  the  center  node.  The  model  is  now  set  up  so  that  the  following  variables  can 
easily  be  changed  in  order  to  find  the  optimal  solution: 

•  Piezo  size 

•  Piezo  placement 

•  Electric  field 

•  Waveform 

•  Amplitude 

•  Piezo  material 

•  Piezo  coupling 

•  Control  scheme 


2.  Conclusion  and  Future  Work 


The  option  was  initially  created  to  build  and  test  a  Kevlar  composite  plate  with  piezos  integrated 
into  the  plate.  The  tests  evaluated  the  optimal  solution  as  determined  by  the  modeling  of  Phase  I. 
As  stated  earlier,  problems  with  obtaining  reliable  data  reduced  the  effectiveness  of  the  modeling 
in  Phase  I.  Some  of  the  option  work  was  completed  by  testing  integrated  piezos,  but  again,  the 
impact  of  the  tests  was  reduced  by  the  unavailability  of  reliable  displacement  data. 

Efforts  are  under  way  to  improve  the  displacement  data.  The  models  have  been  created  and  are 
ready  to  be  optimized.  The  period  between  the  end  of  Phase  I  and  the  start  of  the  option  will  be 
used  to  improve  the  displacement  data.  This  will  enable  Mide  to  complete  the  optimization  very 
early  in  the  option,  so  that  the  prototypes  can  be  built  and  tested  without  disrupting  the  program 
further. 
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Appendix  A.  MATLAB  Main  Code 


%  Construct  a  lumped  paramater  model  of  a  bullet  impact 
close  all 
clear  all 
pack 

%  SIMULATION 
Delt_Frac=0.02; 

%  Delt_Frac=0.1; 

Nstep=4000; 

zeta=2e-10;  %  Damping  ratio 

%  zeta=0*le-17;  %  Damping  ratio 

iplt=l;  %  Plot  option  (0  =  No  plots) 

number_node=0;  %  =  1:  Plot  the  node  numbers 
number_element=0;  %  =  1:  Plot  the  element  numbers 
Nskip=20;  %  Number  of  steps  between  animation  plots 

Deformation_Scale=l;  %  Scaling  factor  for  the  deformation 
Case_Name='APG'; 

save_opt=l;  %  =  1;  Save  the  plots  in  jpg  format 

movie_opt=l;  %  =  1;  Save  a  movie  of  the  impact  response 

stress_opt=l;  %  Calculate  the  stresses 

%  GEOMETRY 

Nx  =  13;  %  Number  of  elements  in  the  planar  direction 

Ny  =  Nx; 

Lx  =  0.2286;  %  Lenght  or  width  of  the  armor  protection 

Nz  =  5;  %  Number  of  elements  in  the  thickness  direction 

Lz  =  .00635;  %  Thickness 

Ly  =  Lx;  %  Consider  a  square  plate 

Ly  =  Lx;  %  Just  a  strip  to  better  match  the  stiffness 

MidN  od=f  ix(Nx/2) ; 

IimpadHVlidNod+MidNod^Nx+1; 

Ibehind=Iimpact  +  Nx*Ny*(Nz-l); 

delx=Lx/(Nx-l); 
dely=Ly/ (Ny-1 ); 
delz=Lz/(Nz-l); 

%  BULLET 

Bullet_Velocity=600;  %  Bullet  impact  velocity 
Bullet_Mass=.00156;  %  Bullet  mass 

%  MATERIAL  PROPERTIES 
Ey  oungXX=24 . 1 E9; 


41 


Ey  oung  YY =24 . 1 E9; 

Ey  oungZZ=l  0 .4E9; 
nhuzx=0.12; 
nhuyz=0.12; 
nhuxy=0.4; 

Density=1440; 

fprintf('\n  Total  Mass  =  %6.2f  kg  or  %6.2f  lb  \n',Density*Lx*Ly*Lz,Density*Lx*Ly*Lz*2.2); 

%  Construct  the  coordinates 
ii=0; 

for  k=l:Nz 
for  j=l:Ny 
for  i=l:Nx 
ii=ii+l; 

Coord(ii,l)=(i-l)>f'delx; 

Coord(ii,2)=(j-l)*dely; 

Coord(ii,3)=(k-l)*delz; 

end 

end 

end 

Nnode=ii; 

%  Construct  the  "Nodal"  Connectivity 
%  First  in  the  X-direction 
ielem  =  0; 
for  k=l:Nz-l 
for  j=l:Ny-l 
for  i=l:Nx-l 

ielem  =  ielem+1; 

Nconec(ielem,l)=(k-l)*Nx*Ny  +  (j-l)*Nx  +  i; 

Nconec(ielem,2)=(k-l)*Nx*Ny  +  (j-l)*Nx  +  i  +  1; 

Nconec(ielem,3)=(k-l)*Nx*Ny  +  j*Nx  +  i  +  1; 

Nconec(ielem,4)=(k-l)*Nx*Ny  +  j*Nx  +  i; 

Nconec(ielem,5)=k*Nx*Ny  +  (j-l)*Nx  +  i; 

Nconec(ielem,6)=k*Nx*Ny  +  (j-l)*Nx  +  i  +  1; 

Nconec(ielem,7)=k*Nx*Ny  +  j*Nx  +  i  +  1; 

Nconec(ielem,8)=k*Nx*Ny  +  j*Nx  +  i; 
end 
end 
end 

Nelem=ielem; 

Ali  ve=ones  (Nelem,  1 ); 

%  Plot  the  grid  on  a  non-equal  axes  grid 
if  (iplt~=0)  ;plot_model_8;end 
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NGP  =  3; 

%  Now  calculate  and  construct  the  finite  element  model 
KMat=sparse(zeros(3*Nnode,3*Nnode)); 

TMat=sparse(zeros(3*Nnode,3*Nnode)); 

MMat=sparse(zeros(3*Nnode,l)); 
for  ielem=l:Nelem 
for  j=l:8 

XI(j,l)  =  Coord(Nconec(ielem,j),l); 

YI(j,l)  =  Coord(Nconec(ielem,j),2); 

ZI(j,l)  =  Coord(Nconec(ielem,j),3); 

end 

[STIFM1  ,JAC]  =  ANSTFSOL  (XI  ,YI  ,ZI  ,EyoungXX  ,EyoungYY  ,EyoungZZ 
nhuxy  ,nhuyz  ,nhuzx  ,NGP  ); 

[EMASS  ,JAC  ,XCG  ,YCG  ,ZCG  ,TMASS]  =  MASSOL  (XI  ,YI  ,ZI  ,Density  ,NGP  ); 


for  i=l:8 

IC  =  Nconec(ielem,i); 
for  kl=l:3 
IC2=(IC-l)*3+kl; 

ICS  =  (i-l)*3  +  kl; 
for  j=l:8 

IR  =  Nconec(ielem,j); 
for  k2=l:3 

IR2  =  (IR  - 1  )*3  +  k2; 

IRS  =  (j  - 1)*3  +  k2; 

KMat(IR2,IC2)  =  KMat(IR2,IC2)  +  Alive(ielem)^STIFMl  (IRS  ,ICS  ); 
end 
end 
end 
end 


for  i=l:8 

ID  =  Nconec(ielem,i); 
for  kl=l:3 
ID2=(ID-l)*3+kl; 

MMat(ID2)  =  MMat(ID2)  +  EMASS  (i ); 
end 
end 
end 

for  i=l:Nnode*3 
TMat  (i, :  )=KMat(i, : )  /MMat(i) ; 
end 

%  TMat(l:3,l:3) 

[phi,eval]=eig(full(TMat)); 
[Freq,iord]=sort(sqrt(real(diag(eval)))/(2>{'pi)); 
for  i=l:Nnode*3 
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phil(:,i)=phi(:,iord(i)); 

end 

phi=phil; 
clear  phil 

delt=Del  t_F  r  ac/max  (F  r  eq) ; 

fprintf('\n\n  Lowest  Natural  Frequency:  %7.3f  Hz  \n\n',Freq(7)); 

save  Restart  TMat  KMat  MMat  phi  Freq 

wold=zeros(Nnode*3, 1 ) ; 
wvold=zeros(Nnode*3, 1 ) ; 

KE_Bullet=0.5*Bullet_Velocity*Bullet_Velocity*Bullet_Mass; 
Imparted_Velocity=sqrt(KE_Bullet*2/(MMat(Iimpact*3)  +  Bullet_Mass)); 

wold(Iimpact*3)Hmparted_Velocity*delt; 

wvold(Iimpact*3)=0; 

MMat(Iimpact*3)HVlMat(Iimpact*3)  +  Bullet_Mass; 
MMat(Iimpact*3-l)HVlMat(Iimpact*3-l)  +  Bullet_Mass; 
MMat(Iimpact*3-2)=MMat(Iimpact*3-2)  +  Bullet_Mass; 
for  i=l:Nnode*3 
TMat  (i, :  )=KMat  (i, :  )/MMat  (i) ; 
end 

wvel=(wold-wvold)/(delt) ; 

%  KE0=0.5*wvel' *MMat'*wvel 

if(iplt~=0);plot_deformed_8;end 

w_store=zeros(Nstep,Nnode*3); 

N  stress=N  s  tep/N  skip ; 

VMStress_store=zeros(Nstress,Nelem); 

wvel_store=w_store; 

istress=0; 

for  itime=l:Nstep 

Time(itime)=(itime-l)*delt; 
wvel=(  wold- wvold)/ (delt) ; 

wnew=+2*wold  -  wvold  -  delt*delt*TMat*wold  -  delt*delt*zeta*TMat*wvel; 

w_stor  e  (i  time, : )- wne w ; 
wvel_store(itime,:)=wvel; 
wvold=wold; 
wold=wnew; 

%  Calculate  the  Kinetic  Energy 

KE  (itime)=0 .5*' w vel ' .  *MMat'  *  w  vel ; 

PE(itime)H)T25*(wold+wvold)'*KMat*(wold+wvold); 
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%  Calculate  the  stresses  around  the  impact  point 
if(mod(itime,Nskip)  ==  0); 

fprintff  Time  Step  %d  \n',itime); 
if  (stress_opt  ==  1) 
istress  =  istress  +  1; 
for  ielem=l:Nelem 
for  j=l:8 

XI(j,l)  =  Coord(Nconec(ielem,j),l); 

YI(j,l)  =  Coord(Nconec(ielem,j),2); 

ZI(j,l)  =  Coord(Nconec(ielem,j),3); 

110*3-2,1)  =  wnew(Nconec(ielem,j)*3-2); 

110*3-1,1)  =  wnew(Nconec(ielem,j)*3-l); 

U(j*3,l)  =  wnew(Nconec(ielem,j)*3); 
end 

[SIGX  ,SIGY  ,SIGZ  ,SIGXY  ,SIGYZ  ,SIGZX  ,SIGVM  ,JAC]  =  ... 

ANSTRSOL  (XI  ,YI  ,ZI  ,U  ,EyoungXX  ,EyoungYY  ,EyoungZZ  ,nhuxy  ,nhuyz  ,nhuzx  ,NGP  ); 
VMStress_store(istress,ielem)=SIGVM(9)/lE6;  %  In  MPa 
%  if(ielem==l); 

%  U(l:16)' 

%  VMStress_store(istress,:) 

%  pause 

%  end 

%  if(ielem==4); 

%  U(l:16)' 

%  VMStress_store(istress,:) 

%  pause 

%  end 

end 
end 
end 
end 

save  Restart  TMat  KMat  MMat  phi  Freq  KE  PE  wvel_store  w_store  Time  VMStress_store 

figure(2) 

subplot(211) 

plot(Time,w_store(:,Iimpact*3));xlabel('Time  (seconds) ');ylabel('Z-Displacement  of  Impact  Point  (m)') 
subplot(212) 

plot(Time,wvel_store(:,Iimpact*3));xlabel('Time  (seconds)');ylabel('Z-Velocity  of  Impact  Point  (m/s)') 
figure(2) 

figure(3) 

subplot(211) 

plot(Time,w_store(:,Ibehind*3));xlabel('Time  (seconds)');ylabel('Z-Displacement  of  Back  (m)') 
subplot(212) 

plot(Time,wvel_store(:,Ibehind*3));xlabel('Time  (seconds) ');ylabel('Z- Velocity  of  Back  (m/s)') 
figure(3) 
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figure(4) 

subplot(311) 

plot(Time,KE);xlabel('Time  (seconds)');ylabel('Kinetic  Energy  (J)') 
subplot(312) 

plot(Time,PE);xlabel('Time  (seconds) ');ylabel('Potential  Energy  (J)') 
subplot(313) 

plot(Time/PE+KE);xlabel('Time  (seconds) ');ylabel('Total  Energy  (J)') 
figure(4) 

if(iplt~=0);animate_filled_8;end 
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VALC ARTIER  ARMAMENTS  DIV 
R  DELAGRAVE 
2459  PIE  XI  BLVD  N 
PO  BOX  8800 

CORCELETTE  QUEBEC  GOA  1RO 
CANADA 

2  DEUTSCH  FRANZOSISCHES 
FORSCHUNGSINSTITUT  ST  LOUIS 
H  ERNST  HLERR 

CEDEX  5  RUE  DU 
GENERAL  CASSAGNOU 
F  68301  SAINT  LOUIS 
FRANCE 

1  DIEHL  GMBH  AND  CO 

M  SCHILDKNECHT 
FISCHBACHSTRASSE  16 
D  90552  ROTBENBACH  AD 
PEGNITZ 
GERMANY 
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1  DYNAMEC  RSCH  AB 
A  PERSSON 

PO  BOX  201 

S  151  23  SODERTALJE 

SWEDEN 

2  ETBS  DSTI 

P  BARNIER  M  SALLES 
ROUTE  DE  GUERAY 
BOITE  POSTALE  712 
18015  BOURGES  CEDEX 
FRANCE 

1  FEDERAL  MINISTRY  OF  DEFENCE 
DIR  OF  EQPT  &  TECH  LAND 
RUV  2 
DHAUG 
POSTFACH  1328 
53003  BONN 
GERMANY 

4  FRANHOFER  INSTITUT  FUR 
KURZZEITDYNAMIK 
ERNST  MACH  INSTITUT 
VHOHLER  E  STRASSBURGER 
RTHAM  KTHOMA 
ECKERSTRASSE  4 
D  79  104  FREIBURG 
GERMANY 

1  MINISTRY  OF  DEFENCE 
DGA/SPART 

C  CANNAVO 

10  PLACE  GEORGES  CLEMENCEAU 
BP19 

F  92211  SAINT  CLOUD  CEDEX 
FRANCE 

2  HIGH  ENERGY  DENSITY  RSCH  CTR 
VFORTOV  GKANEL 
IZHORSKAY  STR  13/19 
MOSCOW  127412 

RUSSIAN  REPUBLIC 

1  INGENIEURBURO  DEISENROTH 
F  DEISENROTH 
AUF  DE  HARDT  33  35 
D  5204  LOHMAR  1 
GERMANY 

1  INST  OF  CHEMICAL  PHYSICS 
S  RAZORENOV 
142432  CHERNOGOLOVKA 
MOSCOW  REGION 
RUSSIAN  REPUBLIC 


NO.  OF 

COPIES  ORGANIZATION 

7  INST  FOR  PROBLEMS  IN  MATERIALS  SCI 
S  FIRSTOV  B  GALANOV  O  GRIGORIEV 

V  KARTUZOV  VKOVTUN  YMILMAN 

V  TREFILOV 

3  KRHYZH ANO V SKY  STR 
252142  KIEV  142 
UKRAINE 

1  INST  FOR  PROBLEMS 
OF  STRENGTH 

G  STEPANOV 

TIMIRY  AZEVSKAYA  STR  2 

252014  KIEV 

UKRAINE 

3  INST  OF  MECH  ENGR  PROBLEMS 

V  BULATOV  D  INDEITSEV 

V  MESCHERYAKOV 
BOLSHOY  61  VO 

ST  PETERSBURG  199178 
RUSSIAN  REPUBLIC 

2  IOFFE  PHYSICO  TECH  INST 

E  DROBYSHEVSKI  A  KOZHUSHKO 
ST  PETERSBURG  194021 
RUSSIAN  REPUBLIC 

1  K&W  THUN 
WLANZ 

ALLMEND S STRAS SE  86 
CH  3602  THUN 
SWITZERLAND 

1  R  OGORKIEWICZ 
18  TEMPLE  SHEEN 
LONDON  SW  14  7RP 
UNITED  KINGDOM 

1  MAX  PLANCK  INSTITUT  FUR 
EISENFORSCHUNG  GMBH 

C  DERDER 

MAX  PLANCK  STRASSE  1 
40237  DUSSELDORF 
GERMANY 

2  NATL  DEFENCE  HDQRTRS 
PMO  MRCV  MAJ  PACEY 
PMO  LAV  A  HODAK 
OTTOWA  ONTARIO  KIA  OK2 
CANADA 
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1  OTO  BREDA 
M  GUALCO 
VIA  VALDIOCCHI  15 
I  19136  LA  SPEZIA 
ITALY 

5  RAPHAEL  BALLISTICS  CTR 
M  MAYSELESS  Y  PARTOM 
G  ROSENBERG  Z  ROSENBERG 

Y  YESHURUN 
BOX  2250 
HAIFA  31021 
ISRAEL 

1  ROYAL  MILITARY  ACADEMY 
E  CELENS 

RENAISSANCE  AVE  30 
B  1040  BRUSSELS 
BELGIUM 

1  ROYAL  NETHERLANDS  ARMY 
J  HOENEVELD 

V  D  BURCHLAAN  3 1 
PO  BOX  90822 

2509  LS  THE  HAGUE 
NETHERLANDS 

1  DEFENCE  MATERIEL  ADMIN 
WEAPONS  DIRECTORATE 

A  BERG 

S  11588  STOCKHOLM 
SWEDEN 

2  SWEDISH  DEFENCE  RSCH  ESTAB 
DIVISION  OF  MATERIALS 

SJ  SAVAGE  J  ERIKS  ON 
S  172  90  STOCKHOLM 
SWEDEN 

3  SWEDISH  DEFENCE  RSCH  ESTAB 
L  HOLMBERG  B  JANZON 

P  LUNDBERG 
BOX  551 
S  147  25  TUMBA 
SWEDEN 

1  TECHNION  INST  OF  TECH 

FACULTY  OF  MECH  ENGINEERING 

S  BODNER 

TECHNION  CITY 

HAIFA  32000 

ISRAEL 


NO.  OF 

COPIES  ORGANIZATION 

3  TECHNISCHE  UNIVERSITAT 
CHEMNITZ  ZWICKAU 
ASCHROEDTER  L  KRUEGER 
L  MEYER 
POSTFACH 
D  09107  CHEMINITZ 
GERMANY 

2  TNO  PRINS  MAURITS  LAB 
HPESKES  RIJSSELSTEIN 
LANGE  KLEIWEG  137 
PO  BOX  45 
2280  AA  RIJSWIJK 
THE  NETHERLANDS 

6  CENTRE  DE  RECHERCHES 
ET  D’ ETUDES  D’ARCUEIL 
D  BOUVART  C  COTTENNOT 
SJONNEAUX  HORSINI 
S  SERROR  F  TARDIVAL 
16  BIS  AVENUE  PRIEUR  DE 
LA  COTE  D’OR 
F  941 14  ARCUEIL  CEDEX 
FRANCE 

1  CDR  EUROPEAN  RSCH  OFFICE 

USARDSG  (UK) 

S  SAMPATH 
PSC  802  BOX  15 
FPO  AE  09499-1500 
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